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Abstract 

We discuss analogies between the charge- and spin response functions of 
itinerant magnets. We show that the spin-analog of the charge stiffness is 
not given by the usual spin stiffness ps, but by the dynamic spin stiffness 
Dg, which is obtained from the dynamic spin conductance Gs{uj) in the limit 
of vanishing frequency uj. The low-frequency behavior of Gs{uj) is used to 
define ideal spin conductors, normal spin conductors, and spin insulators. 
Assuming diffusive spin dynamics, we show that the spin diffusion coefficient 
is proportional to lim^^^^o R-6Gs(lij). We exploit this fact to develop a new 
extrapolation scheme for the spin diffusion coefficient in the paramagnetic 
phase of the Hubbard model. 
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Typeset using REVT^ 



I. INTRODUCTION 



In a classic paper on the Theory of the Insulating State W. Kohn pointed out that 
the behavior of the charge stiffness D,, can be used to distinguish the conducting from the 
insulating state of interacting electrons. Physically Dc can be identified with the weight 
of the Drude peak of the frequency- dependent conductivity a{uj). A system with a finite 
value of Dc deserves to be called an ideal conductor, because then a{ijj) diverges as Dd {ioj) 
for a; — »• 0. In contrast, a normal conductor does not have a Drude peak {Dc = 0), and 
the conductivity o^{uj) approaches a non-zero value for vanishing frequency. Finally, the 
insulating state can be characterized by i^c = and a{uj = 0) = 0. Note that Dc is defined 
in terms of the conductivity, a dynamic quantity. Nevertheless, Dc can be obtained without 
explicitly calculating a{uj): as shown by Kohn the second derivative of the free energy 
with respect to a fictitious vector potential (which is equivalent to a twist in the boundary 
conditions on the many-body wave-function) is proportional to the charge stiffness. More 
recently, the arguments of Kohn |l| have been sharpened by Shastry and Sutherland , who 
calculated Dc exactly for Heisenberg-Ising and Hubbard rings using the ansatz of Bethe. 

Another important observable, which is closely related to Dc, is the superfiuid stiffness 
Pc (the motivation for our slightly unconventional notation will become obvious below). 
Whereas Dc characterizes the normal conducting properties of the system, a finite value 
of Pc indicates long-range superconducting correlations. Dc and pc can both be obtained 
as different limits of the wave-vector and frequency-dependent current response function 
Kai3{(i,uj), which is given by the Kubo formula P]. We shall briefly summarize the relevant 
definitions in Sec.|I|. 

While Dc and pc in Hubbard and related models have recently been studied numerically 
by several groups the analogous quantities Dg and ps that characterize the spin dy- 

namics have not received much attention. In fact, some authors seem not to be fully aware 
of the physically different meaning of the spin analog Ds of the charge stiffness on the one 
hand, and the spin analog Ps of the superfiuid stiffness on the other hand. The quantity 
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Ps is commonly called the spin stiffness^ and is related to the change in energy due to a 
small static twist in the directions of the spins at the boundaries of the system 0. A finite 
value of Ps in the thermodynamic limit is a manifestation of (quasi)-long-range magnetic 
order. Obviously, the spin stiffness ps and the superfiuid stiffness pc are analogous quanti- 
ties: both measure the degree of off-diagonal long-range order in the system. But what is 
the spin analog of the charge stiffness, and what is its physical meaning? In this paper we 
shall not only answer this question, but also discuss the more general concept of the spin 
conductance Gs{uj). We show that the limiting behavior of Gs{oj) for vanishing frequency 
uj can be used to define ideal spin conductors, normal spin conductors, and spin insulators, 
in complete analogy with the charge transport. Normal spin conductors are of particular 
interest, because systems with diffusive spin dynamics fall into this category. In Sec.|^ we 
shall use the concept of the spin conductance to develop a new extrapolation scheme for 
calculating the spin diffusion coefficient of these systems. For localized spin models (such 
as the Heisenberg model) we have discussed the spin conductance and its usefulness for the 
calculation of the spin diffusion coefficient in Ref. [0. 

For definiteness we shall consider the repulsive Hubbard model, although our consid- 
erations are easily generalized for other itinerant magnets. The Hamiltonian is given by 
H = f + V, with 

d 



T=-tEE[4cr+a.+h.C.] , (1) 
r a=l 

V = UY.\cUc.^-]^\cUc.l-\] , (2) 
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where r labels the sites of a (i-dimensional hyper-cubic lattice, and are the primitive 
vectors with length a in direction a = 1,. . . ,d. The operators cj^, a =t, i create spin-cr 
electrons at lattice site r, and cj = [Cr|,4j are two-component operators. 
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II. CHARGE RESPONSE: CHARGE AND SUPERFLUID STIFFNESS 



Let us first recall the Kubo formula for the conductivity, which measures the linear charge 
response to an electromagnetic field. The usual frequency-dependent conductivity can be 
written as aa/siuj) = j^a?~''-Gap{uj), where the dimensionless conductance |Q Gctp{i^) can be 
obtained from the current response function /^^^(q, cu), 

Kap{ci,uj + zO+) 



lim 



The current response function has two contributions, 



(3) 



(4) 



where the diamagnetic part D'^^^ is proportional to the thermal expectation value of the 
kinetic energy operator. 



D'^''' = —{-f) = —ypJn\-f\n) . 



dN 



(5) 



and the paramagnetic contribution is for general complex frequency z given by 



-i 1 ,11 lyy-i 



+ 



{n\Ja{q)\m){m\J^{-q)\n) 
{n\Jp{-q)\m){m\Ja{q)\n) 



(6) 



Em — En + hz 

Here \n) denotes a complete set of exact eigenstates of H, andp„ = {Y.m e~^'^^'^)~~^e~^"^'^ are 
the thermal occupation probabilities of states with energies En, where T is the temperature 
(measured in units of energy). The current operators are 

t 



</a(q) 



2i 



-«qr 



h.c. 



(7) 



The charge stiffness tensor [-De] a/3 and the superfiuid stiffness tensor [pc]af3 can be defined 
via the following limiting procedures. 



[Pc]al3 = - 



q-*0 



lim KadiQ, oj + iO^ 

q^O 

hm Kasi^,, ijJ + ^0^ 



(8) 
(9) 



For completeness, let us also introduce the q-dependent superfluid stiffness 0, 

[Pc(q)]a/3 = - lini ir„/3(q,u; + iO+) , (10) 



which probes the response to time- independent transverse electromagnetic fields (see below). 
Eqs.(|^), (P|)-(pi]|) relate physical quantities characterizing the charge dynamics to the appro- 
priate limits of the linear response function Kai3{(i,uj). For convenience, we have chosen a 
gauge where the scalar potential is set equal to zero, so that the electric field is represented 
by a vector potential, E(t) = —c^^dA{t)/dt. This gauge, which is used very often in the 
literature, has the advantage that the current response can be expressed entirely in terms of 
the current- current correlation function. Of course, the physical current response is gauge 
invariant, see Ref. |jlO . 



Quite generally, at long wavelengths (|q|a ^ 1) the current response tensor Ka/3{(i,uj) 
can be decomposed into a longitudinal part {K\\) and a transverse part {K±), i.e. 

Kc,f3iq,uj) = qaqf3Kii{q,uj) + {6^/3 - qaqi3)K±{q,uj) , (11) 

where = ■ q, with = a^/a and q = q/g. The corresponding decomposition for the 
tensor [pc(ci)]a/3 contains only a transverse component, 

[Pc{q)]ap = {Sap - qaqp)pc{q) ■ (12) 

The longitudinal part vanishes. Physically, this is due to the fact that a static longitu- 
dinal vector potential cannot induce any current If we set q = in Eq.(0) we have 



K±{0, u) = K\\{0, u), because for a spatially uniform field the longitudinal and transverse re- 
sponse are identical. Thus, the conductance tensor Gap^u) in Eq.(^ and the charge stiffness 
tensor [-De] a/3 in Eq.(^ are proportional to the unit matrix (for a system with cubic symme- 
try). The usual charge stiffness and superfluid stiffness pc can be identified with eigen- 
values of the corresponding ternsors (|,|]), i.e. [Dc\af3 = Sa/sDc, and [pc\ai3 = {Sa/3 - qaqf3)Pc- 
Note that with our normalization both quantities have units of energy, p^ is proportional to 
the density of the superconducting electrons. The finite value of pc in a superconductor is 



closely related to the screening of the magnetic field, i.e. the Meissner effect. Note that in 
a normal metal pc = 0, which is a consequence of the fact that in the normal metallic state 
static magnetic fields are not screened. Finally, let us point out that the physical meaning 
of the different order of limits in Eqs.(^ and is easy to understand from the Maxwell 
equation cq x E(q, cj) = co'B(q, cu): If we first take the limit u ^ 0, the electric field E 
vanishes while the magnetic field B can remain finite - in this way we probe the Meissner 
effect. On the other hand, if we first let q ^ 0, we are left with an electric field. 



III. SPIN RESPONSE: DYNAMIC AND STATIC SPIN STIFFNESS 



The generalization of the above definitions for the spin degrees of freedom is straightfor- 
ward. Denoting by cr*, i = x,y,z, the Pauli matrices, in d dimensions we may define 3 x d 
spin-current operators, 

t 



a 



x,y,z 



- Cr-a„) - h.c. 

a = 1, . . . , (i 



(13) 



In complete analogy with Eq.(^ we define the retarded spin-current response function 



(14) 



where the diamagnetic contribution V^^^ is given in Eq.(|^), and the paramagnetic term 
P^-^(q, 2;) is simply obtained from Eq.(|) by replacing — >■ and Jp — > Jg. As discussed 
by Chandra, Coleman, and Larkin [|12|, the spin-current response function /('^^^(q, tu) gives 
the spin response to a fictitious vector potential 5y4^(q, cj), which describes a space- and 
time-dependent modulation in the local spin-density. The proper definition of the dynamic- 
and static spin stiffness tensor is now evident, 



[D 1'^' 



af3 



lim 



limi^^^^(q,w + «0^ 

q^O 



(15) 
(16) 
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Furthermore, in analogy with the q-dependent superfluid stiffness introduced in Eq.(|TOD, let 
us define the q-dependent spin stiffness tensor [|T4| 

[Ps{o,)]% = -Y^mK%{c,u; + iQ^). (17) 

Finally, in analogy with Eq.(|^) we introduce the dimensionless spin conductance 

[G.MU-hm— ^^-^-^ • (18) 

Because the operators J^(0) define the ferromagnetic spin-currents, the limit q — in 
Eqs.([l5|), (|l^) and ([TsD implies that we are looking at ferromagnetic correlations. In the 
case of antiferromagnetism we simply should consider the limit q ^ 11 instead, where 11 = 
[vr/a, . . . , vr/a] is the antiferromagnetic ordering wave- vector. A summary of the analogous 
quantities characterizing the charge- and the spin dynamics is given in Table |. We would 
like to emphasize that the spin analog of the Drude weight Dc is given by the dynamic spin 
stiffness Dg, and not by the static spin stiffness ps- It seems that the dynamic spin stiffness 
has not been discussed in the literature on the Hubbard model. Following the terminology 
used for the charge dynamics, a system with Dg > can be called an ideal spin conductor. 
For Dg = and Gs{0) ^ the system is a normal spin conductor, and the spin insulator 
can be characterized by Ds = Gs{0) = 0. 

IV. SPIN DIFFUSION 

A. Spin diffusion coefficient and Thouless number 

To see the connection between the spin conductance and spin diffusion, consider the 
dynamic structure factor for the spin degrees of freedom, 

S'^iq, u;) = ^ E Pn^iEm -En- huj) 

^ n,m 

X {n\S'\m){m\SL^\n) . (19) 
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Here the Fourier components of the spin operators are S'q = J2r^ cJ^Cr- General hy- 
drodynamic arguments |jl3| show that the diffusive spin dynamics manifests itself in the 



following long-wavelength and low-energy form of the dynamic structure factor, 

S-'M = 2^ %_%,r,.^$^.y . (20) 

where we have assumed cubic symmetry and spin-rotational invariance. Here V is the spin 
diffusion coefficient, and the static susceptibility x is 

X = ^^^,r ^S^'i^^^) = ^T.{SoSl) . (21) 

On the other hand, the real part of the dimensionless spin conductance (|18|) is given by 

Re[Gs{u)]% = -n6{hu;)[D,]%+[G:{uj)r^, , (22) 



where the weight of the 5-function can be identified with the dynamic spin stiffness (|15|), 
and the paramagnetic contribution is 



[G,[u)\^p - hm 



1 _ I 
= TT hni — V pJ{Em -En- fiuj) 

x(n|j;(q)|m)(m|J;^(-q)|n) . (23) 

The matrix elements of the current operator = [J^, J^, J^] can be related to the matrix 
elements of the spin operators Sq via the Heisenberg equation of motion. Using the fact 
that the Hubbard interaction is spin-rotationally invariant, it is easy so show that to leading 
order in q ■ a^, 

^h^ = [S^,T] = Y.i^.a^)3M ■ (24) 

Hence, 

d 

{Em - En){n\Sq\m) = ^^(q " aa)(n|J«(q)|m) . (25) 



a=l 



Substituting this expression into Eq.(|T9D, it is easy to show that 



KHC = i^|^^lim(^S''(q,.). (26) 



Assuming now the diffusive form pO| ) of tlie dynamic structure factor, we obtain for tlie 
spin diffusion coefficient 

Moreover, it is not difficult to sliow tliat in tlie presence of spin diffusion the dynamic 
spin stiffness Dg vanishes due to a perfect cancellation between the dia- and paramagnetic 
contributions in Eg. (pISD . Thus, the existence of spin diffusion means that the system is a 
normal spin conductor. 

Eq.(p7|) can be rewritten in a form which emphasizes a deep connection between spin 
diffusion and charge diffusion in disordered electronic systems |T5| , p!6| |7[] . Defining the rescaled 
dimensionless spin conductance p] 



and the energies 



limK(o;)]L , (28) 

UJ — >\j 



E.. = jj , A. = — , (29) 



(where L = aN^^"^ is the linear size of the system) we obtain from Eg. (p^) 

9s = ^- (30) 

This expression should be compared with the well-known Thouless formula g = -Exh/^ for 
the dimensionless average conductance of a disordered electronic system. Here the so-called 
Thouless energy i?Th is defined as in Eg. (^9]) (with V now given by the average charge 
diffusion coefficient of the disordered system), and A is the average level spacing at the 
Fermi energy. Thus, Egs.(p7D and ( pOD are nothing but the Thouless formula for the spin 
diffusion problem. The dimensionless number Qs defined in Eg.(|30D is the corresponding 
Thouless number. In analogy with disordered electrons, a system with Qs ^ 1 can be called 
a spin metal. 
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B. Spin diffusion in the Hubbard model 



The above analogies are not only interesting from a formal point of view, but also useful 
in practice. We now show that Eg. (p?!) offers a new and physically transparent extrapolation 
scheme for directly calculating the spin diffusion coefficient of the Hubbard model. See Ref. 
1^ for a similar calculation for the Heisenberg model, and Ref. 117] for an alternative method 



to calculate the spin diffusion coefficient in the two-dimensional t — J-model. 
After some straightforward manipulations, Eq.(P7D can be cast into the form 

^ = ^ dsCis), (31) 



Tx Jo 

where the correlation function C{s) is given by 



C{s) = ^{[i{s) + i{-s)]i). (32) 



Here I{s) = e^^^Ie and the (dimensionless) current operator / is 



2i V 



J2 sm{k^a) [4^Ckt - 4|Ck J , (33) 



k 



where Ck = N~^^'^ J2r e~^^ ^Cj.. The bracket in Eq.(P^) denotes thermal average with respect 
to the interacting Hamiltonian H = T + V, see Eqs.(|ip and (0). Because the kinetic energy 
operator T commutes with the current operator /, for U = the integral in Eq.(PTD does 
not exist. Then our model is an ideal spin conductor. This is not surprising, because the 
diffusive dynamics in a system without disorder must be a correlation effect. We would like 



to emphasize that Eq.(pID has been derived under the assumption that the spin dynamics 
is diffusive. The divergence of the integral for f/ = simply indicates that in this case our 
assumption is not correct. 

Because for U = the integral in Eq.(^) is infinite, we expect that for small U the 
spin diffusion coefficient diverges with some power of t/U. Of course, for finite U the 
correlator C{s) cannot be calculated exactly, so that we have to make some physically 
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motivated approximation. A standard approximation scheme, which has proven to be quite 
rehable for the calculation of the spin diffusion coefficient of the Heisenberg model at high 
temperatures ||^JT8[1, is the Gaussian extrapolation of the short-time expansion of C{s) to 
long times. Expanding C{s) in powers of s, 



=0 (2^)! 

the coefficients can be written in terms of multiple commutators. Because C{—s 
C{s), only even powers of s appear. The first two coefficients are 



(34) 



C2 = {I 



H 



(35) 



(36) 



where we have used [I,T] = and (/) = 0. Assuming that the higher coefficients are 
consistent with a Gaussian, the long-time extrapolation is 



C{s) ^ Coexp 



Then we obtain from Eg. (pl|) 



hV 



Tx 2 



(37) 



(38) 



Note that so far we have not assumed that the interaction is small, so that Eqs. (p5D-(p8D 
are valid for arbitrary U. For simplicity, let us now evaluate the coefficients Cq and C2 in 
the limit U ^ 0. Then the averages in Eqs.(^) and (^) are easily calculated with the help 
of the Wick-theorem. Specializing to the case of half filling, we obtain after a lengthy but 
straightforward calculation 



Co = 2A,(r) , 

C2 = 4U'A,{T)Bd{T) 



(39) 
(40) 



where 
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MT) = ^ E sin\Ka)f{e^/T) [1 - /(e^/T)] , (41) 

B^iT) = ^ E /(^k/T) [1 - /(ek/T)] . (42) 

Here f{x) = [e^ + 1]~^ is the Fermi function, and the non-interacting energy dispersion 
in d dimensions is et = 2(it7k, with 7k = Y.'i=i cos(k ■ aa). Away from half filhng we 
should replace Ck ^ Ck — /i, where /i is the chemical potential. In the non-interacting limit 
Tx = \Bd{T), so that we finally obtain from Eq.(p8D in the limit U <^t 



% . 2^.^; . (43) 



Note that V diverges for f/ ^ 0, in agreement with the fact that without correlations there 
is no spin diffusion. Obviously, the spin diffusion coefficient cannot be calculated by naive 
perturbation theory in powers of U. 

In the limit T ^ t we may use Ad{oo) = \ and Bd{oo) = |, so that Eg. (^31) reduces to 

hv ^t^ , , 

— = 2^- , T = oo , 44 
U 

independent of the dimensionality of the system. Recall that this result has been derived 
in the weak coupling limit. More generally, at T = oo it is easy to see from Eqs.(PT|), (P^, 
and (0) that hT>/a^ is proportional to t^/U for all values of U. This follows from the 
fact that the expansion (|34D of C{s) is actually an expansion in powers of {UsY, because 
the current operator / commutes with the kinetic energy operator T. Assuming that the 
integral is convergent, we may scale out the [/-dependence by re-defining the integration 
variable s' = Us. This leads trivially to the energy scale t'^/U. We would like to point 
out that at T = oo the spin diffusion coefficients of the spin 5=1/2 quantum Heisenberg 
antiferromagnet and the half filled Hubbard model at strong coupling are not identical, 
although both are proportional to t'^/U (see Refs. [pJ|,|7|). The reason is that for T = oo 
the value of V in the Hubbard model is determined by states with energies larger than 
U, while the mapping to the Heisenberg model is only valid for energy scales smaller than 
U. Only in the interval /U <^ T <^ U the half filled Hubbard model at strong coupling 
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can be expected to have the same spin diffusion coefficient as the corresponding Heisenberg 
antiferromagnet with exchange couphng J = At'^/U. 

Let us now discuss the low-temperature hmit T <^t. Using the fact that for T ^ 

/(ek/T) [1 - /(ek/T)] ^ ^5(7k) , (45) 



it is easy to see that Eq.(|^) reduces to 

'2dt 



hV ^ <- ad 
—r = 2\/tt 



a2 ^ 53/2 



T 

where the numerical constants and are 



1/2 ^2 



- , r«t, (46) 



= T7Hsin^(^^")^(7k) , (47) 
k 

&<i = ^E%k) , d^2. (48) 
k 



In d = 2 the integral in Eq.(^8D is logarithmically divergent (for 00), so that at low 

temperatures 62 is given by 

62 ~ 4 ln(4t/T) . (49) 

Because we have assumed that the system is in the paramagnetic state, Eg . (^61) should 
be valid for temperatures above the magnetic ordering temperature T/y. Keeping in mind 
that in d < 2 there is no long range order at any finite temperature, and that in d > 2 
the ordering temperature Tjy is exponentially small at weak coupling, Eq. (|46|) describes the 
low-temperature behavior of the spin diffusion coefficient in a wide range of temperatures 
that are small compared with the band-width 4dt. Although the precise numerical value 
of the prefactor in Eqs.(^) and ( ^6D depends on our Gaussian extrapolation scheme, the 
energy scale t'^/U in Eq.(^3]) and the low-temperature behavior given in Eq.(PB|) should be 
independent of the details of the extrapolation method. 



V. SUMMARY AND CONCLUSIONS 



In this work we have used analogies between charge and spin response functions of 
itinerant magnets to clarify the concept of the static and dynamic spin stiffness. Starting 
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from the general Kubo formula for the relevant linear response functions, we have shown that 
the dynamic spin stiffness Dg is the precise spin analog of the charge stiffness Dc- The usual 
(static) spin stiffness ps is the spin analog of the superfiuid stiffness pc, and is only finite 
in the presence of long-range magnetic order. Considering the fact that the charge stiffness 
has recently received a lot of attention it is rather surprising that the corresponding 

quantity Ds has not been studied. Partially, this might be due to the fact that Ds has often 
been confused with the static spin stiffness ps- Table | summarizes analogous quantities. 
We would like to emphasize again that a finite value of Dg does not imply the existence 
of long-range magnetic order. A simple example is the Hubbard model for U = 0, where 
Dg > but Ps = 0. A value Dg > simply means that the system is an ideal spin conductor, 
so that the spin transport is not diffusive. The analogy with charge transport is obvious: an 
ideal conductor has a finite charge stiffness > 0, implying an infinite conductivity. But 
a perfect conductor is not necessarily a superconductor. Only in the latter case pc > 0. 

The low-frequency behavior of the dynamic spin conductance Gs{uj) can be used to clas- 
sify the spin dynamics into three categories: ideal spin conductors, normal spin conductors, 
and spin insulators. In Sec.^ we have further analyzed a particular class of normal spin 
conductors, namely systems with diffusive spin dynamics. In this case the dynamic spin 
stiffness vanishes, but the spin conductance Gs{uj) has a finite limit for vanishing frequency, 
which is proportional to the spin diffusion coefficient. The concept of the spin conductance 
and the associated Thouless number offers a new and physically transparent extrapolation 
scheme for calculating the spin diffusion coefficient. 

We hope that our work will stimulate further research in this field. Numerical calculations 
of Gs{uj) and Ds in strongly correlated itinerant magnets are called for. In particular, by 
varying some external parameter (such as temperature, doping, or interaction strength), it 
might be possible to observe transitions between the three types of spin transport discussed 
above. Numerical calculations of the temperature-dependence of the spin conductance might 
also give evidence for spin-charge separation in strongly correlated systems: Very recently 
Si ||20| pointed out that the existence of spin-charge separation manifests itself in different 
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temperature-dependencies of the spin- and charge conductances. 
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TABLES 



TABLE I. Analogous quantities characterizing charge- and spin dynamics. 

Charge Spin Physical meaning 

charge response function _ft'Q^(q,a;) spin response function K'J^{ci,lu) linear current response to an external 

vector potential 

conductance Gaijiy^) spin conductance [Gs(a;)]^^ response to a time-dependent, spa- 

tially constant vector potential 

q-dependent superfluid stiffness^ [Pc(q)]a/3 q-dependent spin stiffness [ps(q)]^^ response to a time-independent, spa- 
tially varying vector potential 

superfluid stiffness pc spin stiffness ps probes long-range correlations (super- 

conducting or magnetic) 

charge stiffness'' Dc dynamic spin stiffness Ds finite for ideal (charge or spin) 

conductor 

^In the book by Tinkham P| the transverse eigenvalue of this tensor is denoted by K{(\). 
^Also known as Drude weight 
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